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Abstract 

We study plasma effects on radiative transitions (e.g., decay of excited states of atoms or atomic 
nuclei) in a dense plasma at the transition frequencies uj ^ ujp (where ujp is the electron plasma 
frequency) . The decay goes through four channels - the emission of real transverse and longitudinal 
plasmons as well as the emission of virtual transverse and longitudinal plasmons with subsequent 
absorption of such plasmons by the plasma. The emission of real plasmons dies out at uj < ujp, but 
the processes with virtual plasmons strongly enhance the radiative decay. Applications of these 
results to radiative processes in white dwarf cores and neutron star envelopes are discussed. 
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I. INTRODUCTION 



Radiative processes in stars are very important. First of all, they determine heat transport 
in radiative zones of the stars [l], as well as the radiative transfer and structure of stellar 
atmospheres together with the formation of spectra of stellar radiation . In ordinary 
stars (at the main sequence or around) typical radiation frequencies are much higher than 
the electron plasma frequency oOp of stellar matter. As a result, plasma effects do not affect 
strongly radiative processes. 

However, in a dense matter of compact stars (white dwarfs and neutron stars) the plasma 
frequency can be higher or comparable to characteristic radiation transition frequencies u, 
and the plasma effects cannot be ignored. For instance, in a strongly degenerate nonrela- 
tivistic electron gas at a density of 10'^ g cm~^ (that is typical for degenerate cores of white 
dwarfs and outer envelopes of neutron stars) one has hup ^ 0.6 keV. In this example the 
plasma effects can easily affect radiative transitions in atoms and ions. In the inner crust of 
a neutron star at a density of 10^^ g cm~^, where the degenerate electrons are ultrarelativis- 
tic, the plasma frequency becomes very large, hup ~ 3 MeV (depending on the composition 
of the crust; e.g., Ref. also see Sec. IVTl) . This is large enough to influence radiative 
transitions in atomic nuclei. 

The importance of plasma effects for radiative processes in dense stellar matter has been 
mentioned in the literature (e.g., Ref. M). In particular, the plasma effects on the radiative 
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thermal conductivity have been studied in Refs. [6|, 17| but these studies are not fully complete 
(Sec. IVip . For another example, consider an emitter (an ion or atomic nucleus) in an excited 
state with the transition frequency u to the ground state that is lower than Up. What will 
happen with this emitter taking into account that radiative transitions with the emission of 
any electromagnetic quanta are now forbidden? Will it live at the excited level forever? 

These questions can be answered using the available theory of electromagnetic transitions 
in a plasma. The plasma impact on electromagnetic transitions in an non-relativistic labo- 
ratory plasma and a rarefied non-relativistic cosmic plasma has been studied for a long time 
(e.g., Refs. j^, 9]). The plasma effects can modify the emission of electromagnetic quanta 
|9[]. Moreover, collective plasma processes open another electromagnetic transition channel - 
the emission of virtual plasmons and successive absorbtion of these plasmons by the plasma 
jsl, [lo|. The most pronounced of these effects seems to be collisional broadening of energy 
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levels of atoms and ions and associated broadening of spectral lines. It can be important in 
cosmic and laboratory plasmas [10]. 

To the best of our knowledge, the theory of radiative transitions in a plasma (e.g., Refs. 
[sl, [lo|) has been correctly applied only to study non-degenerate laboratory and cosmic plas- 
mas. In this paper we investigate the radiative transitions in a dense degenerate relativistic 
electron gas, particularly at u < Up. The paper is organized as follows. In Sec. [TTlwe outline 
the formalism for calculating electromagnetic transitions rates in a dense plasma. It is sim- 
ilar to the formalism of stopping power for a charged particle moving in a plasma In 
Sec. Uniwe outline the main properties of a degenerate electron plasma. Section [IV] is devoted 
to the radiative decay in the plasma for those transitions that are allowed in the electric 
dipole approximation. In Sec. |V]we address similar problem for the electric quadrupole and 
magnetic dipole transitions. In Sec. IVII we discuss the main results and some applications, 
particularly, for calculating radiative thermal conductivity in white dwarf cores and neutron 
star envelopes and for studying radiative decay of excited states of atomic nuclei and kinetics 
of neutrons in the neutron star crust. We conclude in Sec. IVIII 



II. FOUR RADIATIVE TRANSITION CHANNELS 

Let us consider an external emitter (for instance, an atom or atomic nucleus) immersed 
in a plasma. The plasma is assumed to be uniform and isotropic; it is characterized by the 
longitudinal and transverse dielectric functions ei{u!,k) and etr{uj,k), respectively. We are 
interested in the transition rate Wi^f [s^^] at zero temperature (T = 0) from an upper state 
i to a lower state / whose energy separation is huj. The expression for Wi^f can be written 
as m 
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= J J ^^"^ 

(47r) 



e" ^ /■ /■ \jfi{k)-kf \jf^{k)xkf 



+ 



(1) 



uj'^e\{uj, k) uj'^etv{oi!, k) — k'^c^ 

Any elementary transition (characterized by the given energy loss hu) is accompanied by 
the transfer of an elementary excitation with a wavevector to a plasma coupled to elec- 
tromagnetic field. The integration is performed over all allowed values of k (with k = \k\, 
and dfifc being a solid angle element in the direction of k). Furthermore, 

jf,{k) = J dVjfi{r)exp{-tkr) (2) 
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TABLE I: Four transition channels in plasma environment. 



Channel Plasmon Open at Comment 

Atr real transverse ujp < u Dominates at u ^ uj, 

Al real longitudinal < uj — uip ^ ujp 
Btr virtual transverse any iv 

Bl virtual longitudinal any uj Dominates at uj ^ lj. 



is the Fourier transform of the local transition current jfi{r) [12]. The latter current can be 
calculated from relativistic theory with stationary (relativistic) wave functions of the emitter 
in the i and / states. The transition rate is cumulative. It includes contributions 
of transition channels with different k and different structures of plasma-electromagnetic 
field excitations. Particularly, it intrinsically contains a sum over polarizations of emitted 
plasmons (see below). It neglects the contribution of two- and multiple-plasmon processes 
which is expected to be small in the cumulative rate. 

The dielectric functions (wit h sp atial dispersion) in Eq. ([1]) take into account the plasma 



id ]. In particular, Eq. ([T]) includes the contribution of direct 



effects on the transition rate 
interaction of the emitter with plasma electrons. 

The transition rate ([T]) can be decomposed as (Table [T]) 

w.^f = wUf + wt^f = wfif + wfif + wt^f + wftY (3) 

Here, wj^j and w^L^f correspond to the longitudinal and transverse channels [the terms in 
([T]) containing ei{u!, k) and Etriuj, k), respectively]. Each of these terms, in turn, contains two 
contributions - (A) the emission of a real longitudinal (Al) or transverse (Atr) plasmon and 
(B) the emission and absorption of a virtual longitudinal (Bl) or transverse (Btr) plasmon. 
Not all of the four channels can be opened at once (see below). 

The emission of real plasmons (channels Al and Atr) is allowed in the presence of the 
poles in the denominators of Eq. ([1]), that is at 

ei{uj, k) = 0, u'^EtriuJ, k) = k'^c^. (4) 

The roots of these equations give the plasmon dispersion relations k\{uj) and kx_r{oj) for 
longitudinal and transverse plasmons, respectively. The emission rates for real longitudinal 



and transverse plasmons are then given by the standard expressions 
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J dVlk \jft{k) • k\ 



(47r) 



Atr 



e 



j dVtu \3ft{k) X k\ 



(47r) 



k=ki(ijj) 

d [a;^etr(t^, k) — k'^c^] 



dk 



-1 



fc = fctr{(^) 



(5) 
(6) 



In principle, there can be several poles for one uj; then one should sum over the poles in 
these equations. 

The processes Bl and Btr in Eq. ([1]) involve virtual plasmons. These processes are allowed 
if the dielectric functions have imaginary parts for some values of A; at a given u. Nonva- 
nishing imaginary parts of dielectric functions ensure that the plasma can directly absorb 
electromagnetic fluctuations induced by the emitter. From Eq. ([1]) one has [l^ 
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Bl 
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dk I dfifc 



w 



Btr 
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(47r) 



\jf^{k) ■ k 


^ 5j£:i(co', k) 




ei{uj,k)\^ 



dk / dfifc 



(47r) 



\jfi{k) X kf'^etr{uj,k) 
|a;^£tr('^, k) — k'^c^f 



(7) 
(8) 



Virtual plasmons do not obey any specific dispersion relation and can have a wide spectrum 
of wavenumbers k for a given u. Equations ([7]) and ([H]) describe the effects which are 
in common with collisional broadening of spectral lines (and associated enhancement of 
radiative transition rates) in atomic physics 10|. 

Equations ©-([H]) can be used to study radiative transition rates of a relativistic emitter 
in any uniform and isotropic dispersive medium (we set T — » and disregard thus induced 
transitions). 

In vacuum, where £i = = 1, only one channel survives out of the four. It is Atr - 
the emission of real transverse plasmons, and these plasmons become identical to ordinary 
photons. Then Eq. ([6]) reduces to the well known expression 
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7,,Atr _ vac 



dQi 



jfi{k) X k 



(9) 



(47r) 



where k = k/k. 

Let us stress that the existence of the four radiative decay channels and general expres- 



sions for the partial decay rates have been known long ago (e.g., Refs 



). However, 
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this formalism has been mostly applied to radiative transitions in a non-degenerate plasma, 
where radiation frequencies are typically much higher than Up. In this case the exchange 
of virtual plasmons is usually unimportant, and the authors focused on the emission of real 
plasmons that was not strongly affected by the plasma environment. We will apply the above 
formalism to analyze dense degenerate electron stellar matter where the plasma effects are 
pronounced much stronger. 



III. PLASMA ENVIRONMENT OF DEGENERATE ELECTRONS 



Let us study plasma effects in a strongly degenerate (zero temperature) ideal electron 
gas of any degree of relativity in the absence of a magnetic field. We will comment on the 
effects of finite temperature, ion plasma polarization and magnetic fields in Sec. IVII We 
employ the collisionless dielectric functions etrl"^, k) and e\{uj, k) of a relativistic electron gas 
at T = derived by Jancovici [13] in the random phase approximation. We do not present 
his cumbersome expressions here (note that they should be corrected [l^ at certain values 
of oj and k) but discuss their main properties relevant to our study. 

The most important quantity is the electron plasma frequency. 



iUp = ^Ane'^ne/ml, (10) 
where is the electron number density, m* = yUe/c^ is the effective electron mass on the 



Fermi surface, /ig = y^lC^ + c'^p'f the electron chemical potential (electron rest-mass energy 
included), Pf = h{3n'^ney^^ being the electron Fermi momentum. We also introduce the 
electron Fermi velocity vp = PF/ml- 

The functions etr{i^,k) and ei{uj,k) are generally complex. Their real parts describe 
plasma effects on the propagation of electromagnetic fluctuations, while their imaginary 
parts describe dissipation of such fluctuations. Under typical parameters in dense stellar 
matter for the processes of our study (Sec. [T]), the main source of dissipation is provided 



by the Cherenkov-type absorption at uj < kvp (e.g., Ref. ll|). At T = the dissipation 
switches on abruptly in this domain; it is absent whenever uj > kvp. Thus the integration 
over k in Eqs. ([7]) and ([8]) can be truncated at = uj/vp. Furthermore, in dense stellar 
environment it is reasonable to assume that radiative transition energies are not too large, 
hoj <^ vpPp, and hup <^ vppp. This smallness of hou and hup with respect to typical electron 
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energies greatly simplifies the consideration. 

At u; ~ cijp and k <^ ujjvp the dissipation effect is absent, and the dielectric functions 
take the form 



,2 /, ,2 



e.KA:) ^ 1-^1 + ^^ , (11) 



£tr(^,fc) ^ l-O^^K- (12) 

From Eqs. (jl]) and (|T2l) we immediately obtain the dispersion relation for the transverse 
waves 

= + c^k\ (13) 

This is a good approximation for all k. It corresponds to two transverse plasma modes with 
different polarizations, but the same dispersion relation. The wave frequency satisfies the 
inequality Utj. > kvp at any k. Therefore, these waves undergo no collisionless damping. 
One has c^tr — cUp as ^ 0; and tutr ~ fee as A; ojp/c. In the latter case these waves turn 
into ordinary photons which are almost unaffected by the plasma environment. 

From Eqs. (jl]) and ffTTj) one can derive the dispersion relation for the longitudinal (electron 



Langmuir) plasma waves 



■l = ul + --vlk\ (14) 
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This equation is valid at A; ^ Up/vp-, when uj\ is only slightly higher than cUp, and the 
collisionless damping is absent. At higher k the dispersion equation must be solved numer- 
ically. The solution shows that at some uj\ (~ uOp) the derivative de\{uj^ k)/dk becomes very 

,,A1 



large. This means that the transition rate switches off when the transition frequency 



uj exceeds some value (a few Up). 

It is important that the frequencies of longitudinal and transverse plasma waves are 
always higher than ujp. This implies that corresponding transition rates undergo the plasma 
frequency cutoff, 

wf^f = w^if = at to < Up. (15) 

Finally, let us outline absorption properties of degenerate electron plasma. For typical 

I 

conditions in dense stellar matter, there are two domains [l3| in the {lu, /c)-plane, where the 
imaginary parts of the longitudinal and transverse dielectric functions are non-zero. The 
first domain is given by the inequahty hu < fi — Epp^^k at Tik < 2pF] the second domain 



is determined by \Epp_fik — jj.] < fiuj < Ep^j^f^k — jJ' at any fc, with Ep = Jc^p^ + mlc^. The 



analytic expressions for the imaginary parts of the dielectric functions in these two domains 
are different. Although we have used exact expressions in computations, we notice that in 
the ultrarelativistic gas it is sufficient to consider the only one domain ui/vp ^ k < 2pp/h^ 
where, to a good approximation. 



1 



- 1 + 



(16) 
(17) 



with X = kvp/u. 



IV. ELECTRIC DIPOLE TRANSITIONS 

Let us calculate the transition rate Wi^f in the electric dipole approximation (El). The 
approximation is valid at ka <^ 1, where a is a typical size of the emitter, and k a typical 
plasmon wavenumber. In this case, the transition current is independent of k, jfi{k) ^ 
jfi{0). The angular integration gives 



J dn^UM-kf = -^\jf,{o)fe, (18) 

(477) 

J d(]fc \jf,{0) xkf = ^ \3dO)fe. (19) 

(47r) 

Then the in-vacuum transition rate (associated with the emission of ordinary photons) is 



given by the standard expression 



12| 



where we have used the relation j/j(0) ~ —iurji, Vfi being the position- vector matrix 
element. 

It is convenient to rewrite Eq. ([3]) as 

W,^f = wJ^fR = (^Al + RbI + RAtr + RBtr) , (21) 

where Rai, Rbi, Rait, and i?Btr are the factors, which describe the plasma effects on the 
transition rates in the four channels (Table [T]), and R is the cumulative factor. The partial 
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factors are given by 

Rai = ^ Ji{k)e 



c 



R 



Bl 



R 



Btr 



( 

2c^uo 

71 



dei{iu, k) 



Jtv{k)k 



dk 

k=ki{u)) 

d [a;^£tr(i^, k) — k'^c^] 



dk 



dk 



dk 

Ji{k)eQei{uj,k) 
\ei{uj,k)f 

Jtrik)k^^etrii^,k) 
|a;^£tr('^, k) — k'^c^[ 



-1 



k = kty:(uj) 



(22) 
(23) 
(24) 
(25) 



The functions J\{k) and Jti{k) describe non-dipole corrections to the El approximation at 
large k (at /ca ^ 1), 
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Mk) = ^ IjfiiO) 

Mk) = ^ UMo) 



d^k jfiik) • k 



-2 



dQi 



jfiik) X k 



(26) 
(27) 



For ka ^ 1, we have Ji(fc) — > 1 and Jtr{k) — > 1. 

Equations (122|) - (!25|) determine the plasma corrections to the El transition rate. Let us 
calculate them in a degenerate electron gas. 

The emission of real longitudinal and transverse plasmons occurs at k ^ uj/vp- For the 
transitions at frequencies uo not larger than several ujp in a dense degenerate electron gas, 
one typically has k '^pp/h. This means, that one can use the classical dielectric functions 
ei{u!,k) and etr{uj,k) [l3| for calculating Rai and -RAtr- We have checked, that the use of 
the exact (quantum) dielectric functions has no noticeable effect on the results. Similarly, 
because for the emission of real plasmons we typically have A; -C 1/a, we can always neglect 
the non-dipole corrections in Eqs. (1221) and fl2^ and set Ji{k) = Jtr{k) = 1. Because no 
plasmon emission can occur at < Up, the transition rates w^j- and wf^j^ are suppressed 
as — s> Up. Indeed, at u Up the longitudinal and transverse dielectric functions are given 
by Eqs. ^ and Using then Eqs. ([22]) and ([23D, we find 



-Rai 

-^Atr 



1 
2 



1 1 atO<ti; — ci;p<^ cOp, 



UJ. 



UJ 



^ at a; > ujp. 



(28) 
(29) 
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Equation (l29l) remains a good approximation at all uj > Up] in the limit of u ^ Up the 
factor Rait tends to its in- vacuum value, -RAtr = 1. In contrast, Eq. fl28l) is valid only for u 
close to Up. For higher u, the factor _Rai is strongly suppressed by the fc-derivative of the 
longitudinal dielectric function in Eq. fl221) (no longitudinal plasmons can propagate at high 
frequencies. Sec. IIIip . 

The transition rate at u ^ Up arises from the virtual-plasmon channels B, Eqs. (!24l) and 
fl25|) . First of all, consider the factor Rbi, Eq. ([21]). Substituting the Eq. into ([21]), we 
obtain 

R.. =u-y i^y 7 - (^-'44-] (30) 

2\vfJ \ujJ { x\ x1,^) \e,{uj,xu/vf)V 

where Xmax = 2ppVF/{fiuj). The integrand has a logarithmic singularity which is avoided 
owing to a natural integration cutoff at x = Xmax- Therefore, the integral has the meaning 
of a Coulomb logarithm. According to Eq. (I30l) . we can neglect the non-dipole corrections 
and set J\{xuj/vf) = 1 provided 2pF ~ ^/o- If so, -Rbi depends only on plasma character- 
istics (and does not depend on the transition properties of the emitter); in this sense, i?Bi 
becomes universal. In the opposite case of 2pF ^ Tija, the function Ji{xuj/vf) suppresses 
the integrand in Eq. fl30l) at k ^ 1/a. Then the cutoff of the Coulomb logarithm occurs 
at smaller x = VF/{auj) < Xmax, and the universal factor, calculated from Eq. (130!) with 
J\{xujIvf) = 1, gives the upper limit of i?Bi- 

Let us study the universal regime of 2pF ^ h/a and consider the behavior of Rbi at <^ 
Up. In this case we can use the static longitudinal dielectric function in the denominators of 
Eq. or fl30l) . The asymptotic behavior of Rbi is 

Rb\ oc {up/uY at u <^ Up, (31) 

implying a strong enhancement of the transition rate over the in- vacuum rate [although the 
emission of real plasmons is forbidden, Eq. (fT5l) ]. 

The consideration of the factor _RBtr is similar. First of all, substituting Eq. ( fTTl) into 
Eq. ( !25l) we conclude that there is no logarithmic divergency for the transverse channel. 
This is because of the extra k'^c^ term in the denominator of Eq. ( l25i) . An analysis shows 
that the main contribution to i^Btr comes from intermediate values of x (whereas the main 
contribution to -Rbi comes from large x). As a result, non-dipole corrections to -Rbi are 
much less important than to -RBtr? and Jtrixu/vp) = 1 is a much better approximation than 
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Ji{xu!/vf) = 1. The asymptotic behavior of Rbit is 

i?Btr oc (ujp/uj)'^^^ at u; <^ Up. (32) 

Thus, the transitions through the transverse channel are less efficient than those through 

the longitudinal channel. 

Finally, we have set J\ = Jtr = 1 and calculated the total plasma enhancement factor R 

with the precise dielectric functions [l^ . In the case of ultrarelativistic degenerate electrons 

{vp = c) the factor R depends on the only one argument u = uj/uOp. In the interval 

0.01 < M < 20 the numerical results can be fitted by the expression 

^ 3.0316, , , 29.3 (m-1) + (m- 1)3 

0.93 + 35.6 (n — 1) + (m — l)'^ 

where 0(x) is the Heaviside step-function. The maximum fit error ^ 1.6% takes place at 

u = 1.34. In the limit of m -C 1 the plasma factor R is dominated by Rbi', in the opposite 

limit of u ^ 1 it is dominated by -RAtr- 



V. ELECTRIC QUADRUPOLE AND MAGNETIC DIPOLE TRANSITIONS 

Now consider the case in which the electric dipole transition z — / of the emitter is 
forbidden but the electric quadrupole (E2) or magnetic dipole (Ml) transition is allowed. 
Because the El transition is forbidden, we have 

j>(0) = / dVjf.ir) = 0. (34) 

Multipole transitions are given by next terms in the expansion of the transition current over 
ka. For the E2 and Ml transitions, the transition current can be written as 

jfi{k) = J dV j fi{r) exp{-ikr) ^ -i J dV^J/i(^) (^ ' r). (35) 



15| . we obtain 



Using the standard expansion over spherical vectors Ijtv/ 

J=0 M=-J 

where r = r/r. The terms with different J correspond to different transition types. 

The term with J = 2 refers to an E2 transition. Indeed, a rank 2 spherical vector can be 



Jdk) = --- E E / dVrjf.ir) ■ Yj;,{r), (36) 



presented as [15 1 



Y2Mir) = grad {r'YMr)) , (37) 

VlOr ^ ^ 
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where 1^21// (^) is a spherical function. Then one can rearrange the integral term in Eq. (136|) 



as 



j dVrjfi{r) 



lUO 



r 



-1 



lUJ 



57r 



Q 



2-M5 



(38) 



where Q^2-m an electric quadrupole component 



12j and pfi{r) is the matrix element of the 



density operator, that is related to the transition current through the continuity equation 



iujpfi{r) + div j/i(r) = 0. 



The J = 1 term in Eq. (1361) corresponds to an Ml transition. Using the relation 



(39) 



(40) 



one finds 
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J dVrjf.ir) . = ^ / dVrjf.ir) • [r X grad YUr) 



-zf-1 



>M. / 3 „(m) 



(41) 



where Q^^\.j is a component of the magnetic dipole moment. 

Finally, the J = term in Eq. fl5Bl) is non-standard. It is absent in vacuum, but appears 
in plasma. Rearranging the spatial integration in Eq. fl36|) in the same way, as for E2 and 
Ml transitions, we obtain 



ldVrjf,{r).Y,l*ir) 



dVp/i(r) 



Q 



2, 



(42) 



2V47r J ' 2V47r 
where 

Q2 = j dVpf,{r)r\ (43) 

In order to separate the contributions from the above terms to the longitudinal and 
transverse transition channels let us split the fc-dependent spherical vectors into components 
longitudinal and transverse to k. The J = 2 term contains both, transverse and longitudinal, 
components |15[ |: 

Y,\Ak) = ^Y.^lkk) + ^Y,[-,'\k), (44) 

where Y2M^\k) is the longitudinal spherical vector and l^jjj^(fc) is the transverse electric-type 
spherical vector. The J = 1 term in fl36l) contains only the transverse component. 



^lA//(^) — ^Im(^) 



(45) 
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Yi^{r) being the transverse magnetic vector. The J = term in fl36|) contains only the 
longitudinal vector, 



Now we can calculate the transition rate from Eq. (p 
we find for the transverse channel: 



(46) 

Performing angular integration, 



dnk\kx jf^{k)\' 



Stt 



. M 



+ 



20 



2-M 



M 



(47) 



The two terms in the square brackets correspond to the Ml and E2 transitions, respectively. 
No contribution from the J = is present, because the J = term is purely longitudinal. In 
vacuum, only the transverse channel contributes to the transition rate. Then from Eq. ([6]) 
we recover the well-known expression 



Atr 



4a/ 



(m) 



-M 



M 



— "^vac ' vac 



(48) 



The angular integration of the longitudinal part of Eq. ([T]) gives 



4:71 

dQk \k-jfi{k)f = — /cV 



2-M 



L M 



(49) 



The first term in the square brackets corresponds to the E2 transition, while the second term 



refers to a different, purely longitudinal transition [9|. The latter transition is not forbidden 
by the standard selection rules and should be kept in line with the E2 transition (at M = 0). 
Note, that this term can be presented as a trace of the quadrupole moment tensor of the 
emitter, Q2 = TrlQa/j}, if this tensor is defined in a non-standard (not irreducible) form as 
Qai3 = / Pfi{r) XaXp. If SO, the Q2 term can be regarded as an additional contribution to 
the quadrupole transition [q]; nevertheless it can be presented even for a spherical emitter. 
Let us add, that there is no Ml transition in Eq. ( H9l) - it is forbidden because the Ml 
transition current is purely transverse. 

The plasma effects on the transition rate can be described by introducing the plasma 
factors R in accordance with Eq. fl2T|) . We obtain 



E2 



W 



E2 



R 



E2 



W- 



E2 



.Ml 



^vac-f<Ml 



.Ml 



p(2) p(2) p(2)\ 
(-"^Al "I" ^Atr "T ^Bl "I" -"^Btr ) 

p(2) p(2) \ 
^Atr "T -"-Btr ) 5 



(50) 
(51) 
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where 



d(2) 



(2) 
Atr 



d(2) 



d(2) 

^BtT 



TTCiJ 



c}ei(u;, k) 
d [uj'^etT{uJ, k) — k'^c^] 



dk 



dk 

k^^ei{LU, k) 



fc = fctr(l^) 



|2 ' 



dfc 



\e\{uj,k)\' 

k^^St^iu), k) 
lu'^EtriyJ, k) — k'^c 



|2 ' 



(52) 
(53) 
(54) 
(55) 



the upperscript (2) marks the second-order multipole expansion. Equations differ 
from Eqs. fl2^ - fl2^ by powers of k in the numerators {k^ instead of the k"^) and by pre- 
factors. Moreover, the total transition rate in the second-order muhipole expansion contains 
an additional term [9j] in the longitudinal channel, 



where 



w. 



w 



E2 



5\Q2 



4E Q^'^ 



M 



2-M 



(56) 



(57) 



In order to calculate this term one should know transition matrix elements. 

We have calculated and fitted the factors Re2 and Rui under the same assumptions as 
the factor R for the El transitions. We considered an ultrarelativistic degenerate electron 
gas {vp = c) and employed exact dielectric functions [l3[. Again, Re2 and Rui become 
functions of the only one parameter u = uo/ujp which was varied in the range 0.01 <u< 20. 
Our fit to -Re2 is 



R 



E2 



. 607.8,^ 
1 + ^^(1 



0.0048m^ 



(58) 



(2) 

7.65. Note, that the function R^{ does not 



with the maximum fit error of 1.6% at u 
deviate from its small-w asymptotic behavior R^} = 607.8 in the entire fit interval. 
The fit to Rmi is 



= 1 + A (4-82 - 0.71nM - 0.47m) 



(59) 



with the maximum fit error of 3.8% at u = 0.1. 
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FIG. 1: (Color online). Different partial contributions to the total plasma enhancement factor R 
(curve 'tot') as a function of uj/ujp for El transitions at vp = c. Other solid curves show i?Btr) -^Bl) 
i?Ai) and Rait- The dashed curve is the low-w asymptote Rbi = S.OSujp/u)'^. 

VI. DISCUSSION 

In Fig. [T]we plot various plasma factors for an El transition rate as a function of u/up 
in the ultrarelativistic strongly degenerate electron plasma {vp = c) at femaxCt ^ 1- The 
solid line marked 'tot' shows the total plasma enhancement factor R. Other solid curves are 
partial contributions Rbit, Rb\, Rai, and Rait given by Eqs. (I22l) - (p5l) . The factors Rai and 
i?Atr vanish at uj < Up because no real plasma waves can be emitted under such conditions. 
At uj ^ 3ujp the main contribution to the total transition rate comes from the emission of real 
transverse plasmons. In the limit of uj ^ ujp the plasma effects disappear and R ~ -RAtr 1- 
Radiative transitions via virtual longitudinal plasmons always dominate over transitions via 
virtual transverse plasmon, R-qi > -Retr- All transitions aX uj < Up go via the exchange 
of virtual plasmons, the transition rate being greatly enhanced in comparison with its in- 
vacuum value. The dashed line in Fig. [1] shows the low-a; asymptote Rbi = 3.03 (up/uj)'^; it 
is accurate at u < Up, where Rbi dominates. 

All quantities, plotted in Fig. [H are calculated neglecting non-dipole corrections to the 
transition current [by setting Ji{k) = Jtr(fc) = 1 that is valid at fcmaxfl = ^ppa/h <^ 1, see 
Eqs. fl26|) and fl27|) ]. The factors R in Fig. [1] are universal, and depend only on u/up (in 
the limit of vp ^ c). However, the condition /cmaxO ^ 1 can be violated. Such a violation 
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FIG. 2: Factor C versus A^max'^ an illustration of the effect of non-dipole corrections on /^bi in 
Eq. (j60p for a simplified model of radiative deexcitation of atomic nucleus of radius a (see text). 
The inset shows Ji{ka) for this model. 

does not significantly affect Rai, Rait, and Rbh but can change Rbi (Sec. HV|) . The effect of 
non-dipole corrections on i?Bi is demonstrated in Fig. [2l Now the universality is lost and the 
result depends on the specific form of the function J\{k) (determined by the wave functions 
of the emitter). For illustration, we consider a dipole transition from the lowest excited 
state with orbital momentum L = 1 to the ground state {L = 0) in a spherical potential well 
of radius a with infinitely high walls (as a very rough model of El deexcitation of atomic 
nucleus). The appropriate function Ji{ka) is shown in the inset. At u; ^ cjp we still obtain 
the asymptote 

-Rbi = C'(A;maxa)(^^p/c^)^ (60) 

where Clk^g^^a) is now determined by Ji{ka) (and does not depend of ujp/uj). The function 
C(^maxCt) takes into account the non-dipole corrections. It is shown in Fig. [2] for our partic- 
ular model. The increase of fcmax^ reduces C{k^g_^a) with respect to its purely dipole limit 
C(0) ~ 3.03. A reduction by a factor of 2 is achieved at k^^^a ^ 10. 

Thus, non-dipole corrections lower the transition rate, but the rate remains enhanced 
over its in- vacuum level by a factor of (co'p/ct;)^.This is because the expression for i?Bi (at 
u < Up) contains the factor {up/u)'^ which arises from Qe\{uj,k). The integration over k 
in Eq. ( 130|) can be carried out using static dielectric function e\{0,k); the function Ji{k) 
specifies only a numerical prefactor, but does not violate the {up/uj)'^ dependence. 
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FIG. 3: (Color online). Plasma enhancement factors -Re2 and -Rmi versus Lojuj.p at vp = c, together 

(2) (2) (2) (2) 
with the factors -Rgi , -RBtr' -^Ai ' ™^ ^Atr ^o^' separate transition channels. 

Now let us consider the E2 and Ml radiative transitions. Their principal features remain 
the same as for the El transitions. In Fig. [3] we plot the plasma enhancement factors -Re2 
and -Rmi- One can see that the plasma enhancement of E2 and Ml transitions is stronger 
than for El transitions. However, this is true only if the El transition is forbidden. If not, 
the effect of higher-order transitions (particularly, E2 and Ml) is included in the functions 
Ji(/c) and Jtr(^) and has already been discussed above. 

Other curves in Fig. [3] show partial contributions to the plasma enhancement factors 
from different radiative decay channels. The main contribution to E2 transitions comes 

(2^1 f2^ 

from (as for El transitions). Typically higher values of R^l [with respect to -Rbi for 
the El case, see Eq. flM|) ]. result from the appearance of an additional in the numerator 
of Eq. (15^ . This leads to a stronger cUp/o; dependence in the asymptotic behavior of -Rbi'', 
in comparison with i?Bh ('^p/'^)^ instead of (tUp/a;)^. Note that the asymptotic expression 
= 607.8 (cup/cu)^ remains an excellent approximation in the entire range of uj presented 
in Fig. [3l Note also, that we have used the first nonvanishing term in the series expansion 
of the transition current over ka. Therefore, our results for E2 and Ml transitions are valid 
for l-pp <^ hja. In the opposite case, just as for the El transitions, the results will depend 
on the exact form of the local transition current jfi{r) 

We cannot plot the contribution of the additional term, wf_^j- [see Eq. (1571)]. to the total 
transition rate in a similar universal form. If the corresponding moments were equal, then 
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the transition rate in the longitudinal channel in a plasma at cu ^ cUp would be about twice 
larger than the E2 transition rate. For uj ^ ujp the transition rate wf_^j: vanishes. 

While calculating the i?-factors, we have used the dielectric function of the degenerate 
electron gas and have neglected the ion contribution. This approximation is expected to 
be valid for transition frequencies u which are much higher than the ion plasma frequency 
Upi (see Fig. Hj). Because Upi <^ Up, the transition frequencies u ^ Upi, at which the ion 
contribution can be important, are much lower than the electron plasma frequency. If 
necessary, the ion contribution can be studied using similar approach. 

The same plasma effects occur in a magnetized plasma but the magnetic field complicates 
the problem. Because of the anisotropy, introduced into the plasma polarization properties 
by the magnetic field, the plasma waves (plasmons) become of mixed type (neither longi- 
tudinal, nor transverse) and have many branches (for instance, electron cyclotron modes). 
The properties of the plasmon emission (channels Al and Atr) of an atom in a rarefied 

nn 

magnetoactive cosmic plasma were studied, for instance, in Refs. [9|, 116[. The effect of the 
magnetic field on the processes with virtual plasmons seems to be unexplored. 

In our analysis, we have employed zero-temperature approximation but similar effects 
should be pronounced at finite temperatures. Moreover, thermal plasma fluctuations, avail- 
able in this case p/^], can power inverse transitions and excite the emitter. The efficiency of 
inverse transitions depends on temperature and plasma parameters. 

The plasma effects are important for studying a number of phenomena in neutron stars 
and white dwarfs. These effects are outlined below using the energy-density diagram for 
dense stellar matter (Fig. HI left vertical scale). The solid lines marked as Tp, Up, and 
Upi show the density dependence of the electron degeneracy energy ksTp (Tp being the 
electron degeneracy temperature), as well as the electron and ion plasma energies, hup and 
hupi. For simplicity, we employ the model of accreted neutron star crust (Table A3 from 
Ref. jisl). The curves Tp and oOp are rather insensitive to possible variations of nuclear 
composition in an accreting neutron star (and quite close to the curves for a neutron star 
whose crust is composed of the ground-state matter {4]). The curve Upi is more sensitive to 
the composition but is relatively unimportant for our analysis. Note that the neutron drip 
occurs at p ^ 6 X 10^^ g cm~^ in the accreted crust (and at p ~ 4 x 10^^ g cm~'^ in the 
ground-state matter). Typical temperatures in neutron stars and white dwarfs are below 
10^ K [ksT ^ 0.1 MeV). Degenerate electrons become relativistic at p ^ 10^ g cm~^. 
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FIG. 4: Energy-density diagram (left vertical scale, solid curves) for dense stellar matter. The 
curves Tp, ujp and LO^i show, respectively, the Fermi energy of degenerate electrons, electron plasma 
energy TiWp, and ion plasma energy fiuj^i for an accreting neutron star. The short-dashed vertical 
line positions the neutron drip point. The dotted line (right vertical scale) shows A:niax versus p to 
characterize the importance of non-dipole corrections for El transitions. See text for details. 

Plasma effects can affect beta captures in dense stellar matter, for instance, in the crust 
of an accreting neutron star in a binary system with a low-mass companion. Such systems 
manifest themselves as X-ray transients which demonstrate periods of active accretion and 



quiescence 



19| . Observations show that neutron stars in X-ray transients remain warm 



during quiescent periods that is o: 
nuclear transformations 
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ten explained [20| by deep crustal heating associated with 



22l |. particularly, beta captures, in the accreted matter. 
When the accreted matter is gradually compressed by newly accreted material, the density 
in local matter elements goes up increasing the Fermi energy of degenerate electrons. This 
triggers beta captures with the appearance of daughter nuclei in ground or excited states. 
If the daughter nuclei are born in the excited states, they can de-excite through radiative 
transitions 3|; the associated energy release can contribute to the deep crustal heating. 
The nuclear composition of the accreted matter can be very different and contain a wide 



spectrum of nuclides 2^. This means numerous beta captures involving various nuclei at 
the densities up to 10^^ — 10^^ g cm~^ (Fig- HI)- In this case the electron gas is strongly 
degenerate and ultrarelativistic, the electron plasma energy Tilo^ can reach a few MeV and 
become larger than transition energies Tiuj in some nuclei. What will happen with these 
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nuclei? A naive answer would be that they would not decay to lower states because they 
cannot emit any electromagnetic quanta at < ojp. Our results show quite the opposite. 
The plasma environment enhances the decays through the processes B involving virtual 
(mostly longitudinal) plasmons. The dotted line in Fig. H] plots (right vertical scale) the 
values of A^max- As follows from the above discussion (see Fig. [2]), at /cmax^ ~ 10 our plasma 
enhancement factor for El transitions starts to deviate from the universal enhancement (!33ll . 
The dotted line indicates that, for typical radii a of atomic nuclei, we have fcmaxCt ^ 10 at 
any p in Fig. HJ so that the enhancement remains universal. Let us stress, however, that we 
use a very crude model of El transition in Fig. [21 We would advise to check the condition 
for the breaking of universality (1331) in specific situations. 

Another example is provided by the reactions involving neutrons (n) in accreting neutron 



stars 



231 ]. Specifically, we mean the reactions (n,7) and (7,n) (neutron absorption by a 



nucleus with the emission of electromagnetic quantum, and an inverse process). These 
reactions can occur at densities 10^^ — 10^^ g cm~'^ near the neutron drip density in the 
neutron star crust (Fig. Hj). They can accompany deep nuclear burning of accreted matter 
and affect energy release and nuclear transformations in deep crustal heating process as well 
as X-ray superbursts (highly energetic X-ray bursts demonstrated by some accreting neutron 
stars). Again, many nuclei can be involved, and typical energies huo of electromagnetic 
transitions can be lower than hujp. Our results cannot be used directly to study the neutron 
reactions, but they can be modified for that purpose. They demonstrate that the plasma 
effects cannot suppress j23l] the neutron capture reactions (n,7) at uop > uj. Moreover, we can 
expect that even at ^ ujp, but at not very low temperatures, there will be a substantial 
level of fluctuating plasma microfields (associated with virtual plasmons) to power the inverse 
reaction (7,n). 

Finally, the present results can be useful for calculating the radiative thermal conductivity 
in a degenerate electron gas. This is an important problem for outer cores of white dwarfs and 
outer envelopes of neutron stars, where the radiative conduction becomes comparable to the 
electron one (the latter dominates in the deeper, strongly degenerate layers of these objects; 



see, e.g., Ref. 



24|). With increasing density into the degenerate matter, the electron plasma 



frequency becomes comparable to typical radiative transition frequencies {u ~ ksT/h) and 
then exceeds them. Radiative conduction is provided by real electromagnetic waves (not 
virtual excitations), which leads to the plasma cutoff of the radiative thermal conductivity at 
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low temperatures (/c^T ^ hup, Fig. Hj). This cutoff has been mentioned in the astrophysical 



literature (e.g. 



A general physical theory of radiative transfer in dispersive media was 



constructed long ago [25|. Several attempts have been made (e.g., to calculate the 

radiative thermal conductivity in dense stellar matter with account for the plasma effects. 
However, these calculations have neglected the contribution of longitudinal plasmons which is 
expected to be important at ksT < hup, especially in the non-relativistic mildly degenerate 
electron gas (T ~ Tp) where the radiative thermal conductivity can be comparable with the 
electron one. 



VII. CONCLUSIONS 

We have analyzed the radiative transition rate of an emitter (an atom or atomic nucleus) 
immersed in a dense degenerate plasma. Such a transition goes, generally, through four 
channels which involve real and virtual longitudinal and transverse plasmons (Refs. 9|, llO| ; 
Table [T]). The emission of real plasmons is allowed only at radiative transition frequencies lj 
higher than the electron plasma frequency Up. The processes with virtual plasmons operate 
at any u. 

Our main conclusions are: 

1. The cumulative effect of the plasma is to enhance the radiative decay rate over the 
standard radiative decay rate through the emission of photons in vacuum. In the limit 
of ^ cUp the plasma enhancement effect disappears. 

2. The enhancement becomes especially strong at <^ cUp (where real plasmons cannot 
exist at all), being mainly provided by processes with virtual longitudinal plasmons. 

3. The plasma enhancement takes place for electric dipole transitions, and for higher- 
order transitions (such as electric quadrupole and magnetic dipole one); it is more 
pronounced for higher-order transitions. 

4. In a strongly degenerate ultrarelativistic electron plasma the plasma enhancement de- 
pends mainly on the parameter uo/ujp. This dependence is calculated and approximated 
by analytic expressions for El, E2 and M2 transitions. 
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The plasma enhancement effects can strongly modify radiative thermal conduction in 
dense stellar matter, kinetics of atomic nuclei in excited states, emission and absorption of 
neutrons. Such effects can be important in degenerate cores of white dwarfs and envelopes 
of neutron stars but are almost unexplored. 
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